Diffusion with Stochastic Resetting 
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We study simple diffusion where a particie stochastically resets to its initial position at a con- 
stant rate r. A finite resetting rate leads to a nonequilibrium stationary state with non-Gaussian 
fluctuations for the particle position. We also show that the mean time to find a stationary target 
by a diffusive searcher is finite and has a minimum value at an optimal resetting rate r* . Reset- 
ting also alters fundamentally the late time decay of the survival probability of a stationary target 
when there are multiple searchers: while the typical survival probability decays exponentially with 
time, the average decays as a power law with an exponent depending continuously on the density of 
searchers. 
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'Stochastic resetting' is a rather common process in ev- 
eryday life. Consider searching for some target such as, 
for example, a face in a crowd or one's misplaced keys 
at home. A natural tendency is, on having searched un- 
successfully for a while, to return to the starting point 
and recommence the search. In this Letter we explore 
the consequences of such resetting on perhaps the most 
simple and common process in nature, namely, the diffu- 
sion of a single or a multiparticle system. We show that 
a nonzero rate of resetting has a rather rich and dramatic 
effect on the diffusion process. 

The first major effect of resetting shows up in the posi- 
tion distribution of the diffusing particle. In the absence 
of resetting, it has the usual Gaussian distribution whose 
width grows diffusively ~ \fi with time. Upon switch- 
ing on a nonzero resetting rate r to its initial position, 
this time-dependent Gaussian distribution gives way to a 
globally current-carrying nonequilibrium stationary state 
(NESS) with non-Gaussian fluctuations, given in Eq. ([%]>, 
The process of resetting manifestly violates detailed bal- 
ance and thus provides an appealingly simple example of 
a NESS. 

Resetting also has a profound consequence on the first- 
passage properties of a diffusing particle. The study of 
first-passage problems and survival probabilities of diffus- 
ing particles arises in diverse subjects such as in reaction- 
diffusion kinetics, predator-prey dynamicsllj, as well as 
in persistence in nonequilibrium systems [2j . Such prob- 
lems are fundamental to nonequilibrium statistical me- 
chanics as they involve irreversible processes not obey- 
ing detailed balance. Related models are also relevant 
to the study of search strategies in ecology or sampling 
techniques for the characterisation of complex networks. 
For example, intermittent searches involve diffusive mo- 
tion combined with long range movements of the searcher 
and mimic the scan and relocation phases of foraging an- 
imals 
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Figure 1. Schematic space-time trajectory of a one dimen- 
sional Brownian motion that starts at xq and resets stochas- 
tically to its initial position xq at rate r. 



A well-studied problem is the mean time for a sta- 
tionary target at the origin to be absorbed by a single 
diffusing particle (trap) or a team of diffusing traps dis- 
tributed with uniform density. Many significant results, 
such as the fact that the mean time to find the target by 
a single diffusing particle diverges and that the survival 
probabilty of the target decays with time as a stretched 
exponential in the presence of a finite density of diffus- 
ing particles, are well-established when there is no reset- 
ting 

Our results, summarised in this paragraph, show that 
these statistical properties of the search for a stationary 
target by diffusing traps (searchers) are fundamentally 
altered when the searchers reset stochastically to their 
initial positions. In the case of a single diffusing trap, 
the mean time to find the target at the origin, given by 
Eq. ([7]) , becomes finite in the presence of a non-zero re- 
setting rate r and as a function of r, has a minimum at 
a nontrivial value r* given by Eq. ijHJ. Thus there is an 
optimal resetting rate that makes the target search most 
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efficient. The survival probability of the target decays 
exponentially with nontrivial decay rate summarized in 
Eqs. (|12ll3p . For the case of mutiple traps distributed 
with uniform density p and each diffusing independently 
with a diffusion constant D, the effect of resetting is even 
more dramatic. For example, in one dimension and in 
the absence of resetting, both the typical and the aver- 
age survival probability of the target decay with time t 
as a stretched exponential, ~ exp(— Xpy/D t), albeit with 
two different values of the decay constant A @, [l(| • In 
contrast, we find that when one switches on the resetting 
the typical survival probability of the target asymptoti- 
caly decays faster with time, as an exponential, Eq. (|26p. 
On the other hand, the average survival probability de- 
cays more slowly asymptotically, as a power law with 
a nonuniversal exponent that varies continuously with 
density p, Eq. (|2"5|) . This slower algebraic decay can 
be attributed to contributions arising from rare, extreme 
events. 

We now present a derivation of our results. For simplic- 
ity, we focus below on the one-dimensional case, though 
generalization to higher dimensions is straightforward 
and will be mentioned at the end. To begin with we 
consider the simple case of a single particle diffusing in 
one dimension with diffusion constant D and stochasti- 
cally resetting to its initial position at a constant rate 
r (see Fig. [T|). The Master equation for p(x,t\xo), the 
probability that the particle is at x at time t, having 
begun from xq at time 0, simply reads 



dp(x,t\x ) d 2 p(x,t\x ) 



dt 



D- 



dx 2 



-rp(x,t\xo)+rS(x—x ) (1) 



with initial condition p{x, 0) = 8(x— Xq). The second and 
third terms on the rhs and represent a negative probabil- 
ity flux rp(x, t\xo) out of each point x and a correspond- 
ing positive probability flux into xo which sums to r. It 
is easy to show that in the long-time limit the particle 
attains a stationary distribution 
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exp(-a |^ - x \) 



(2) 



where a = y/r/D is an inverse length scale corre- 
sponding to the typical distance diffused by the particle 
between resets. The distribution ([2]) is evidently non- 
Gaussian, with a cusp at x = xq. Since resetting creates 
a source of probabality at xq while probability is lost from 
all x ^ xq, there is circulation of probability even at long 
times, making (|2|) a non-equilibrium stationary state. 

The result ([2]) has a simple and appealing interpreta- 
tion in terms of a renewal process. The time evolution 
of the diffusive particle (see Fig. ([T])) is indeed a re- 
newal process, where each resetting event to xq renews 
the process at a Poisson rate r and between two consec- 
utive renewal events, the particle undergoes free diffu- 
sion. Consider the particle at some fixed time t. If r de- 
notes the time elapsed since the last renewal, the position 



distribution is simply Gaussian, e~( x ~ x °^ / 4£,T /y/AirDr. 
Clearly < r < t is a random variable and for large t, 
t is simply exponentially distributed, Prob(r) = re~ rT . 
Thus, for large t, the stationary position distribution of 
the particle is obtained by averaging over r 
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which precisely yields the result in ©. 

We now turn to the first-passage properties through 
the origin of the diffusing particle in the presence of reset- 
ting. The origin can be thought of as a stationary target 
which is absorbed when the diffusing particle, starting 
initially at xo > at t — 0, hits the origin for the first 
time. The diffusing particle also resets to its initial po- 
sition xo at a constant rate r. To solve the first-passage 
problem, it is advantageous to use the backward Master 
equation approach where one treats the initial position 
as a variable. Hence, we first set the initial position to be 
at x > 0, different from the resetting position xo, solve 
the problem with arbitrary x and xq and eventually set 
x = xq. Let Q{x,t) denote the survival probability or 
persistence of the target up to time t (it also depends 
implicitly on the resetting position xq). This is then also 
the probability that the diffusing particle, starting at x, 
has not hit the origin up to time t. The backward Mas- 
ter equation (where the initial position x is the variable) 
then reads 



dQ{x,t) ^ D fQ(x,t) 



dt 



dx 2 



rQ{x,t) + rQ(x ,t) (4) 



with boundary and initial conditions Q(0,t) = 0, 
Q(x,0) = 1. The second and third terms on the rhs 
correspond to the resetting of the initial position from x 
to xo, which implies a loss of probability from Q(x, t) and 
a gain of probability to Q(xQ,t). The Laplace transform 
q(x,s) = / °° dt e~ st Q(x,t) satisfies 
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d 2 q(x, s) 
dx 2 



- (r + s)q(x,s) = -1 -rq(x ,s) (5) 



which can be exactly solved by noting that the general 
solution is q(x, s) = Ae ax + Be~ ax + (l + rq(xo, s))/(r+s) 
with a = y/(r + s) /D. Now q(x, s) remains finite as 
noo therefore A = 0. The condition g(0, s) = fixes 
the constant B = —(1 + rq(xQ, s))/(r + s) and q(xo, s) is 
determined self-consistently as 



q(x ,s) 



1 — exp(— axo) 
s + r cxp(— axo) 



(6) 



The mean first-passage time, T(xo), may be obtained 
from T(x Q ) = - f °° dtt 9Qi g t °' t] = q{x ,s = 0) which 
yields 
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T(x ) = - (exp(a a;o) - 1) 



(7) 
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where we recall ao = y/r/D. The first point to note 
is that T is finite for r > and diverges as r — > as 
T ~ r' 1 / 2 , which recovers the well-known result that the 
mean time for a diffusive particle to reach the origin (in 
the absence of resetting) is infinite. Also T diverges as 
r — > oo, the explanation being that as the reset rate in- 
creases the diffusing particle has less time between resets 
to reach the origin. 

Now consider T as a function of r for a fixed xq. Since 
T diverges as r — > and r — > oo it is clear that there 
must be a minimum of T with respect to r at an optimal 
value r* . We consider the dimensionless variable z = 
aoXg, the ratio of the distance of the intial site from 
the target to the length diffused between resets, with 
which our results may be simply expressed. Then the 
condition for the minimum, dT/dr = 0, reduces to the 
transcendental equation 



which has a unique non-zero solution z* = 1.59362 

Thus, there is an optimal resetting rate r* = (z*) 2 D/xq 
in one dimension that minimizes the search time to find 
the target at the origin. 

It is difficult to invert the Laplace transform (|6]) ex- 
plicitly for all parameters. However, one can deduce the 
large t asymptotic behaviour for fixed parameters r, xq. 
Generally this is determined by the singularity, in the 
complex s plane, of q(xo, s) with largest real part. From 
([6]) there will be a pole in q(xo, s) when 

s + rc- ax " = (9) 

which may be rewritten as s — r(u — 1) where u satisfies 

u = I - cxp(-u 1/2 z) . (10) 

An addition to the trivial solution u = (which corre- 
sponds to the branch point at s = — r) there is exactly 
one non-trivial solution u of (|10p where < u < 1. 
Thus q(xQ,s) has a simple pole at sq = —r(l — u ) the 
residue of which determines the asymptotic behaviour as 

Q(x ,t) ~ e s °* . (ii) 

zuq + 2u Q — z 

Different limiting values of uq and hence sq can be de- 
duced. For z -C f , ito ~ z 2 and sq ~ — r + rz 2 whereas 
for z S> 1, uq ~ 1 — e -2 and sq ~ — rcT z . These results 
imply that the survival probability decays asymptotically 
as 

Q(x ,t) ~ 2z 2 cxp(-r£(l - z 2 + ...)) for z< 1 (12) 
Q(xo,t) ~ exp(— r£exp(— z)) for z>l. (13) 

One may relate the large z result to the Gumbel distri- 
bution for the extremum of independent random vari- 
ables [111 as follows. In time t the mean number of resets 



to position xq is N = rt. For the target to survive re- 
quires that between each reset the maximum distance to 
the left attained by the particle is less than xq. In a 
given reset interval of length r, the particle undergoes 
free diffusion around xq- Hence, the probability that the 

minimum is bigger than m is erf ( ^^==£J- J [IJ and the pdf 

of the minimum m is Prob(m) = e^ < - m ^ x "^ l ADr jsfnDr 
with m < xq. Next, we average over t drawn from the 
exponential distribution Prob(r) = re~ rT . Thus, the ef- 
fective pdf of the maximum over each reset interval is 
then p(m) = ao 6 """' 10- " 1 ' where m < xq. Taking the 
iV" intervals as statistically independent, the probability 
that the minimum of all N intervals stays above is 
[J p(m)dm] N . Performing the integral and taking large 
z limit yields w exp [— riexp(— z)\. 

We now turn to the survival probability of a stationary 
target at the origin in the presence of many independent 
searchers (diffusive particles). We consider N diffusive 
particles i = 1, . . . , JV, each of which is reset indepen- 
dently to its initial position Xi with rate r. The survival 
probability of the target is given by 

jv 

P s (t) =Y[Q( Xi ,t) (14) 

i=l 

where Q(xi,t) is the survival probability in the single 
searcher problem. 

The initial positions assumed to be inde- 

pendent and each distributed uniformly over the box 
[— L/2, L/2]. Consequently, P s (t) is a random variable. 
Its average is simply Pf v (i) = (P s (t)) x where () x de- 
notes averages over x^s. However, P s (t) for a typi- 
cal initial configuration is not captured by the average. 
The typical P s (t) can be extracted by first averaging 
over the logarithm of P s (t) followed by exponentiating: 
P s typ (t) = exp [(lnP a (t))a;]. One can draw an analogy to 
a disordered system with P s (t) playing the role of parti- 
tion function Z and disorder variables. Thus the 
average and typical behavior correspond respectively to 
the annealed (where one averages the partition function 
Z) and the quenched (where one averages the free energy 
In Z) averages in disordered systems. 

In the annealed case we get 

Pr(t) = {Q(x,t))Z = cxpiVln[l - (1 - Q) x ) (15) 
where 

i r L/2 

(l-Q)x = T / dx[l - Q(x,t)] . (16) 

Letting N, L — > oo but keeping the density of walkers p = 
N/L fixed, and using the symmetry Q(x, t) = Q(—x,t), 
we obtain 

/•OC 

P s av (i) -> exp-2p / dx[l - Q(x,t)] = cxp-2pM(i) . 
Jo 

(17) 
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The Laplace transform M(s) = J °° M(t) e - St dt can be 
determined using ([6]) 



M(s) 



T + s I s + r 
m 

sra 



It may then be inverted [12] to obtain 



l~ 



1/2 



(18) 



(19) 



the initial conditions whose memory is retained through 
resetting. 

Our results can straightforwardly be extended to 
higher dimensions. For example, in d-dimensions, the 
second derivative term in the Master equation ((TJ) is re- 
placed by a d-dimensional Laplacian. The stationary so- 
lution in d-dimensions is easy to find [l~2T | 



with y = rt, where 

My) = 



y , 1 - e- v 
dv 



erf 



(y - v) 



1/2 



-(y-v) 



- v) 



The function fi(y) behaves asymptotically as 

fj,(y) ~ — p- for small y 
~ In y + 7 for large y 



(20) 

(21) 
(22) 



where 7 is Euler's constant. Thus the long-time be- 
haviour (rt 3> 1) of the annealed survival probability 
is a power law with an exponent that varies continuously 
with the density (fl~7|) as 



P s av (i) ~ At- 2 ^ D ^ 1/2 , 



(23) 



where A is a constant that may be determined. The 
short-time behaviour (rt <C 1) is stretched exponential 
decay P s av (i) ~ ex^-^Dt) 1 / 2 /^) as in the case of 
diffusion without resetting 0, [l(| ■ 

In contrast, the typical behavior (the quenched case) 
P s typ (i) = exp [(mP s (f))a;] can be expressed as 



P? p (t) = exp^i^Qi^t))- = CX P 



2p 



L/2 



In the long time limit, we have using pip 



dx lnQ(x,t) 
(24) 



dxo In Q(xo,t) ~ Constant 



—t I da; |s (a; )| 




(25) 

The integral can be done in closed form [l2j and we find 
J °°da;oSo(a;o) = -(£>r) 1 / 2 4(l - In 2). Thus the asymp- 
totic decay of the quenched total survival probability is 
exponential 



P«(t) - exp ( -t P (Dr)^ 2 8(l - In 2) 



(26) 



The correction to the argument of the exponential in (|26[) 
will come from the branch point at s = — r in ^ and is 
expected to give a contribution C^i 1 / 2 ). The fact that the 
average and typical survival probabilities have distinct 
asymptotic behaviours reflects the strong dependence on 



p st (x\x ) = 



(ao) d 
(27r) d /2 



(a \x- x \YK„(a \x - x \) (27) 



where v = l—d/2 and K v is the modified Bessel function. 
Similarly, the mean absorption time for an absorbing ball 
of radius e at the origin is [ijj 



T(x ) = - 
r 



K v (a Q e) 
x \J K„(a \x \) 



1 



(28) 



The survival probability in case of multiple searchers, 
both annealed and quenched, can also be computed in 
higher dimensions [12 1. 

In conclusion, in this Letter we have shown that intro- 
duction of resetting events leads to rather rich effects on 
simple diffusion process. In particular, the effect of reset- 
ting on the efficient search for a stationary target both 
by a single or a team of searchers is rather profound. Our 
results can be extended in several directions, in particu- 
lar one may consider resetting rates r(x) that depend on 
the position of the particle [12j. It would also be inter- 
esting to study the effect of resetting on other available 
search strategies based on nondiffusive processes such as 
Levy flights |3,[i|. 

We thank R. A. Blythe, B. Derrida and S. Janson for 
useful discussions. 
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